SPECIFICATION 
CONTROL FOR AN I.S. MACHINE 



5 The present invention relates to an LS. (individual section) machine and 

more specifically to a control for such a machine. 

BACKGROUND OF THE INVENTION 

10 The first LS. machine was patented in U.S. Patents Nos. 1,843,159, dated 

February 2, 1932, and 1,911,119, dated May 23, 1933. An LS. (individual 
section) machine has a plurality of identical sections. Each section has a frame on 
which are mounted a number of section mechanisms including blank side and 
blow side mold open and close mechanisms, an invert and neck ring mechanism, 

15 a baffle mechanism, a blowhead mechanism, a plunger mechanism and a takeout 
mechanism. Associated with these mechanisms is process air used for cooling, 
for example. Each of the section mechanisms and the process air have to be 
operated at a selected time in the section cycle. 

In the original LS. machine, devices (valves which operated the 

20 mechanisms and the process air, for example) had to be mechanically turned on 
and off each cycle and the timing process was controlled by a 360° timing drum 
which was a cylindrical drum having a number of annular grooves, one for each 
valve, each supporting "on" and "off dogs for tripping a corresponding switch 
associated with a particular valve. The rotation of this mechanical timing drum 

25 through 360° has always been equated to the completion of one control cycle of 
the machine or section and accordingly men skilled in this art have always 
analyzed machine performance in a wrapped cycle, i.e., one that repeatedly cycles 
from 0° to 360°. When electronic timing replaced the mechanical timing drum, 
devices were turned on and off by an electronic sequencer which replicated the 



wrapped 360° control cycle of the mechanical timing drum. An encoder defined 
the angular location of the electronic sequencer, and electronic switches were 
turned on and off at the same angles as was the case with a mechanical timing 
drum. A very significant development that greatly enhanced the power of the 
5 electronic sequencer was the concept of thermodynamic modes (U.S. Patent No. 
3,877,915) wherein groups of these electronic switches were linked so that they 
could be simultaneously adjusted. These machine controllers allow the user to 
electronically adjust the on/off schedule (angle) for the various valves which 
operate the section mechanisms. This conventional approach does not allow an 

10 operator to directly command the machine to achieve desired forming durations 
(e.g. blank contact time, reheat time). It also does not prevent the user from 
setting invalid or even potentially damaging sequences in which the mechanisms 
collide. Only with considerable experience, and process insight can an operator 
effectively adjust the machine timing with the conventional approach and since 

1 5 skill levels vary greatly, the productivity of the machine can vary substantially. 

OBJECT OF THE INVENTION 

It is an object of the present invention to provide an improved control 
20 system for a glass forming machine which will simplify machine operation and 
facilitate the tuning of the machine for higher productivity. 

Other objects and advantages of the present invention will become 
apparent from the following portion of this specification and from the 
accompanying drawings which illustrate a presently preferred embodiment 
25 incorporating the principles of the invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



Figure 1 is a schematic illustration of one section of an I.S. machine which 
can have one or more of such sections, 
5 Figure 2 is the first part of the Network Constraint Diagram for the Blow 

and Blow process; 

Figure 3 is the second part of the Network Constraint Diagram for the 
Blow and Blow process; 

Figure 4 is the third part of the Network Constraint Diagram for the Blow 
10 and Blow process; 

Figure 5 is the fourth part of the Network Constraint Diagram for the Blow 
and Blow process; 

Figure 6 is the fifth part of the Network Constraint Diagram for the Blow 
and Blow process; 

1 5 Figure 7 is the sixth part of the Network Constraint Diagram for the Blow 

and Blow process; 

Figure 8 is the seventh part of the Network Constraint Diagram for the 
Blow and Blow process; 

Figure 9 is the eighth part of the Network Constraint Diagram for the 
20 Blow and Blow process. 

Figure 10 is network model for branch incidence matrix; 

Figure 11 is an event timing chart for a 360° electronic sequencer which 
controls a section of an I.S. machine; 

Figures 12A and 12B are a network diagram for use to unwrap a wrapped 

25 cycle; 

Figure 13 is a block diagram illustrating the creation of a computerized 
model of a mathematical representation of a network constraint diagram 
unwrapped from a wrapped cycle; 
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Figure 14 is a block diagram illustrating the portion of the computerized 
model which converts wrapped event angles into unwrapped event times; 

Figure 15 is a logic diagram illustrating the operation of a control using the 
computerized model to analyze an unwrapped schedule re constraint violations 
5 such as sequence, collision, or duration violation; 

Figure 16 is a logic diagram illustrating the operation of a control using the 
computerized model to analyze an unwrapped schedule to define the duration of 
the thermal forming processes; 

Figure 17 is a logic diagram illustrating the operation of a control using the 
1 0 computerized model to analyze an unwrapped schedule re the optimization of the 
schedule; 

Figure 18 is a logic diagram illustrating the operation of a control using the 

computerized model to define the event angles for a feasible schedule with 

thermal forming process duration "N" inputs; 
1 5 Figure 1 9 is a logic diagram illustrating the operation of a control using the 

computerized model to optimize an unwrapped schedule; 

Figure 20 is a logic diagram illustrating the operation of a control using the 

computerized model to identify, when a schedule is determined to be feasible, any 

active constraint that restricts further improvement; and 
20 Figure 21 is a logic diagram illustrating the operation of a control using the 

computerized model to minimize the wear on the displaceable mechanisms. 

BRIEF DESCRIPTION OF THE PREFERRED EMBODIMENT 

25 An I.S. machine includes a plurality (usually 6, 8, 10, or 12) of sections 

10. Each section has a blank station including a mold opening and closing 
mechanism 12 having opposed mold supports 14 which carry blank mold halves. 
When these mold supports are closed by a suitable displacement mechanism 16 
which can displace the mold support between open (illustrated) and closed 
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positions and which is displaced by a motor 18 such as a servo motor, discrete 
gobs of molten glass can be delivered to the closed blank mold. The open top of 
the blank mold will then be closed by a baffle of a baffle support 22 which is 
displaceable between remote and advanced positions by a motor (such as a servo) 
5 24. If the section is operating in the press and blow mode, the plunger of a 
plunger mechanism 26 is advanced vertically upwardly into the gob to form the 
parison. Cooling air will be supplied to the plunger via a valve VI . If the section 
is operating in the blow and blow mode, the finish is formed by applying settle 
blow air through a valve V2 in the baffle mechanism 22, and the parison is 

10 formed with the application of counterblow air to the plunger via a valve V3, 
while vacuum is applied to the baffle through a valve V4. 

After the parison is formed, the baffle support is retracted, the mold 
supports are retracted and a pair of neck ring holder arms 30 which are rotatively 
supported by an invert mechanism 31, will be rotated 180° by a servomotor drive 

15 32. The blank station also includes a mold opening and closing mechanism 12 
having opposed mold supports 14 which carry blow mold halves. These mold 
supports are displaced between open and closed positions by a suitable 
displacement mechanism 16 which is displaced by a motor 18 such as a servo 
motor. With the parison located at the blow station, the mold supports are closed, 

20 the neck ring arms are opened to release the parison (each arm is displaceable by 
a pneumatic cylinder (not shown) which is operated with a suitable valve V5), 
the invert mechanism returns the neck ring arms to the blank side (the arms close 
prior to arrival) and a blow head support 34 which is displaceable between a 
retracted position and an advanced position where a supported blowhead closes 

25 the blow mold, is displaced to the advanced position by a suitable motor such as a 
servo 36 to blow the parison into the bottle. This final blow is controlled by a 
valve V6. 

When the bottle is formed, the blowhead is retracted, blank molds are 
opened and a takeout mechanism 38 which is driven by a suitable motor 39, such 
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as a servo motor, is displaced to pick up the formed bottle and carry it to a 
location above a deadplate 40 where it is cooled while suspended and then 
deposited onto the deadplate. In addition to the movement of mechanisms and 
devices, process air to mechanisms, moveable or stationary, may also be 
5 controlled. When the blow molds are closed, mold cooling air is turned on to 
cool the formed bottle. 

Each section is controlled by a computer 42 which operates under the 
control of a 360 degree timing drum (programmable sequencer) which defines a 
finite number of angular increments around the drum at which mechanisms, etc., 

10 can be turned on and off each 360 degree rotation. The control knows the time it 
takes for rotating 360 degrees and this time can be fixed or defined as the duration 
between once per cycle pulses such as pulses originating from the feeder of the 
I.S. machine. Each valve is cycled (turned on and off) and each mechanism is 
cycled within the time of one machine cycle by an electronic timing drum 

1 5 (programmable sequencer) which is part of the computer 42. 

In accordance with the present invention a tool is defined first by making 
an unwrapped cycle constraint diagram for an actual I.S. machine configuration 
and then making a mathematical representation of the unwrapped cycle constraint 
diagram that is capable of automated formulation and solution. "Unwrapped" 

20 means the I.S. is a process cycle beginning with the formation of a gob of molten 
glass by severing the gob from a runner of molten glass and ending with the 
removal of a formed bottle from the blow station. This process cycle takes more 
than one 360° machine cycle of the timing drum to complete (normally 2 times 
360° machine cycles). 

25 Figures 2-9 show a possible Network Constraint Diagram for a 

representative blow and blow process for making glass bottles in an I.S. machine. 
The cycle starts with the shear cut represented by time node zl ("z" and V 
denote a time node). Gob Delivery/M13 (a block containing an "M" represents 
an activity that will move between start and end positions with the direction of 
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motion being indicated with arrows) begins at zl and ends with nl77/e26/n6 (a 
vertically oriented equal sign labeled "e" connecting two nodes indicates that the 
two connected nodes occur at the same time). Gob Delivery/Ml 3 motion is 
subdivided into two submotions: 1. Gob In Collision Zone With Baffle/m2 (a 
5 block containing a "m" represents a submotion) which starts at zl/el/n3 and ends 
at n4. Gob Traverses Blank Mold/m3 which starts at n4/e2/n5 and ends at n6. 

Node zl (shear cut) also has another branch Total Process/dl3 which starts 
at zl/e79/nl75 and ends at nl76/e78/n84 (Figure 9). Derived branches are 
identified with ellipses containing "D" and represent thermal process durations 

10 that are defined as a function of the machine events. 

Figure 2 also shows that Plunger To Loading Position/MPl ("P" means 
prior cycle) must be completed at nl3. Node nl3 is the time when the motion 
Plunger To Loading Position/Ml was completed at nl5 during the prior cycle. 
This is indicated by cycle time branch (Figure 6) that connects nl3 to nl5. The 

1 5 plunger includes an independently moveable thimble and at the end of Plunger To 
Loading Position/Ml, both the thimble and the plunger are up. Node nl77, the 
end of Gob Delivery/Ml 3 must be some time (s2) ("s" along side a pair of closely 
adjacent direction arrows represents some time (a sequence constraint) that will 
pass between connected nodes) after n 1 3 . 

20 Figure 2 also shows node n20 which is the time when Baffle Of£7MP15 

was completed in the previous cycle t2. This is indicated by cycle time branch t2 
that is connected to node n22 (Figure 4) which is the time when Baffle Off7M15 
is completed in the subsequent cycle. Node n20 is connected to nl which starts 
Baffle On/M14 some time (s22) after n20, i.e., Baffle On/M14 cannot start until 

25 Baffle Off/15 is completed. Motion branch Baffle On/M14 ends at node n93. 
Baffle motion is broken down into submotions Baffle Moves To Interference 
With Gob/m4 which starts at nl/e27/n7 and ends at n8 and Baffle On 
Completion/m5 which starts at n8/e3/n9 and ends at nl0/e28/n93. Also shown is 
collision branch Baffle Collides With Gob/cl (collision branches are represented 
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by a squiggly line which is identified by "c") connecting node n4 to n8. This 
means that the gob must be at n4 before or no later than the baffle reaches n8 in 
order to be sure that no collision will occur. 

Figure 2 also shows node n40 which is the time Blank Molds Closed/MP9 
5 the last cycle (n40 is connected to node n55 (Figure 6) which is the end of Blank 
Molds Closed/M9 the present cycle with tl indicating a cycle difference). Blank 
Molds Closed/MP9 was complete at n40 which is some time (s21) before the start 
of Gob Traverses Blank/m3 at n5. 

When the gob is fully delivered in the blanks nl77/e24/n26 Blank 

10 Contact/dl (Figure 3) begins and continues until n25/e25/n28 when the Blanks 
Open/M5. Prior to Blank Contact/dl at time n5/e63/nl83 (the time when Gob 
Traverses Blank/m3 begins) a vacuum valve is opened starting the process branch 
Vacuum Assist/pl3 (process branches are identified with ellipses containing "P"). 
Vacuum Assist/pl3 will continue until nl82 when the vacuum valve will be 

15 closed. This means that at the same time that the gob is traversing the blank, 
vacuum is being applied through the neck ring (before completion of the plunger 
moving to the loading position) to help draw the gob into the neck area of the 
blank and into the neck ring. 

At nl2 which is some time (s5) after the gob is delivered (nl77) and some 

20 time (s3) after the baffle is on (nlO), a compressed air valve is opened to start 
Settle Blow/pl which ends at node nl l/e73/n21/e68/nl55 with the closing of the 
compressed air valve. When Settle Blow/pl ends Settle Blow Vent/plO begins 
and ends at nl9 and Neck Ring Contact/d8 begins and ends at nl54/e69/nl 13/ 
with Neck Rings Opening/m21 (Figure 5). This means that at the completion of 

25 settle blow the gob is in contact with the neck ring and heat is being removed 
from the gob until the neck rings are opened. Baffle To Down/M2 (Figure 2) 
begins at n69, which is some time (si) after nil and ends at n35 (closes top of 
blank for counter blow). At n 172 (Figure 3) which is some time slO after nl77 
when the gob is fully loaded into the blank molds and another time s 1 1 after 
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Blank Cooling/pP7 ended at n 173 during the last cycle (til), Blank Cooling/p7 
begins with the opening of a. valve and continues until nl71 when the valve is 
closed. 

At nl56 (Figure 3), which is time s40 after Vacuum Assist/pl3 ends at 
5 nl82 and which is time s7 after nl9 when Settle Blow Vent/plO is over, the 
Plunger Is Displaced To Counter Blow Position/M3 (the thimble is retracted out 
of the glass), a process that ends at n70 and at the same time (nl56/e70/nl58) the 
glass in the area of the finish, which is in complete contact with the molds, will 
reheat (Corkage Reheat/d9) until nl57/e71/nl60 which is some time (s39) after 

10 n70 and which is some time (s3 6) after n3 5 (the end of Baffle To Down/M2. At 
nl60, Counter Blow/pll starts with the opening of a valve and continues until 
time nl59/e80M81/ when a valve opening a vent in the baffle is opened to allow 
the process Counter Blow Vent/p 12 to start. This process ends at n 180. At time 
nl48, which is some time (s38) after nl59, the Plunger (is displaced) To Invert 

1 5 Position/M4 whereat both the thimble and the plunger are down (this takes until 
nl47). 

At nl49/e66/nl51, which is some time (s37) following the end of Counter 
Blow Vent/p 12 at nl80, the following events begin simultaneously: 1. Parison 
Bottom Reheat/d7 which lasts until nl50/e65/n28 and 2. Baffle Off/Ml 5 (Figure 

20 4)which lasts until n22/e30/n33.. Baffle Off7M15 can be divided into two 
submotions; the first being Baffle Off Clears Interference With Invert/ml 1 which 
begins at nl49/e29/n32 (Figure 3) and ends at n31/e7/n34, and the second being 
Baffle Off Completion (past interference)/ml2 which starts at n34 and ends at 
n33. At n28 (Figure 3), which is some time (s8) after nl49, the following events 

25 simultaneously occur: 1 . Blanks Open (blank molds are opened)/M5 which ends 
at n27 (Figure 4) leaving the bottom of the parison on the bottom plate of the 
blank molds; 2. Reheat (parison)/d4 begins at time n28/el5/n29 (Figure 4) and 
continues until n61/el6/n30 (Figure 6) (some time (sl5) after Blow Head 
On/Ml 8 movement is completed at nlOl) when Final Blow/p2 (Figure 7) begins, 
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ending at n63; and Inverted Reheat/d3 which begins at n28/e8/n38 (Figure 3) and 
continues until n37/e9/n39 (Figure 5) which corresponds to the completion of 
Invert/M6 which began at n24. At n36 (Figure 5), some time (s4) following n37, 
reheat will continue with the parison inverted (Parison Invert Recovery/p4) until 
5 nl7. The invert motion is subdivided into a number of submotions. At the 
beginning of invert displacement (n24/e53/nl53) (Figure 4), there is the 
submotion Invert To Baffle Interference/m40 which ends at time nl52/e67/nl25. 
The next submotion is Invert Baffle Interference To Invert Blowhead 
Interference/m32 which ends at time nl24/e52/nl27. The next submotion is 

10 Invert To Takeout Interference 1 From Blowhead Interference/m3 which ends at 
nl26/e60/nl40 when Invert (moves) To Takeout Interference2/m33 begins 
ending at nl39/e61/nl42. The next submotion is Invert (moves) To Takeout 
Interference3/m38 which begins at nl42 and ends at nl41/e54/nl29. Finally 
Invert Completion/m35 (Figure 5) occurs beginning at N129 and ending at 

15 nl28/e55/n39. 

A number of collision branches are identified. Plunger Collides With 
Invert/c2 (Figure 3) when the plunger M4 is not displaced to the invert position 
before the invert moves (time n!47 vs. time n24). Blanks Collide With Invert/c3 
(Figure 4) when the blanks M5 are not displaced to the open position before the 

20 invert moves (time n27 vs. time n24). A number of other collisions are shown: 
Baffle Collides With Invert/c4 when the baffle mil reaches a selected point 
before n24 and Baffle Collides With Invert/cl8 when the baffle mil reaches its 
fully off position before nl52 when the invert has reached the outer limits of its 
interference zone with the baffle. By dividing the interference zone into more 

25 than one zones, the mechanism can get an earlier start. The blow head and the 
invert will collide cl2 if Blow Head Up MP 19 has not occurred (the last cycle t4) 
before the invert has reached the end of the Invert Baffle Interference To Invert 
Blowhead Interference (time n23 vs. time nl24). 
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Also shown is the motion of the takeout: Takeout Through 
Interferencel/mpl3 (Figure 4) which ends at nl43 (the last cycle/t7); Takeout 
Through Interference2/mp24 which ends at n 144 (the last cycle/t8); and Takeout 
Through Interferernce3/mp36 (Figure 5) which ends at nl45 (the last cycle/t9). A 
5 number of collisions are identified: Takeout Collides With Invert/c 1 3 (Figure 4) if 
the invert reaches Interferencel before the takeout does (nl43 vs. nl26). Takeout 
Collides With Invert/c 17 if the invert reaches Interference2 before the takeout 
does (nl44 vs. nl39). Takeout Collides With Invert/cl6 (Figure 5) if the invert 
reaches Interferences before the takeout does (nl41 vs. nl45). At nl79 (Figure 

10 4) which is some time (s34) after n28, Neck Ring Cooling/p9 begins with the 
opening of a valve and continues until nl78, which is some time s35 before n24 
when the invert/M6 begins to move. 

The blow molds, which were open Mp24 (Figure 4) at time nl4 during the 
last cycle tlO, begin to close at time n98/e56/nl46 which is some time (sl7) after 

15 nl4. Closing has a number of submotions: Molds Close to Ware Width/m39 
(Figure 5) which starts at nl46 and ends at nl09/e62/n85; Molds Close To 
Parison Width/ml6 which starts at n85 and ends at n62/e32/n42; Molds Close to 
Receive Position/ml4 which starts at n42 and ends at n41/el0/n44; and Molds 
Close Shut/ml 5 which starts at n44 and ends at n43/e3 l/n97 (Figure 6). Takeout 

20 Clears Ware From Mold/Mp30 (Figure 4) must have operated in the prior cycle t3 
before Molds Close To Ware Width/m39 to avoid a collision of the takeout with 
the molds clO (time n89 vs. time nl09. Furthermore, Parison Invert Recovery/p4 
should be complete before the molds close (time n!7 vs. n62/e32/n42) or a 
collision Parison Collides With Mold/c5 will occur. 

25 The neck rings open to release the parison at the blow head (Neck Rings 

Open/M8) (Figure 5). This motion which occurs from n46 to n45/e44/nll2 is 
divided into two parts: Neck Ring Open Delay/ml 8 which starts at the same time 
n46/e45/nl 1 1 and ends at nl 10/e43/nl 13 (some time (s26) after n41 - the end of 
Molds Close To Receive Position/ml4 and some time (s25) before Blow Molds 
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Close/Ml 6 is complete at n97) (Figure 6) when the second part (Neck Rings 
Opening/m21) starts. This second part ends at nll2. In the event that Neck 
Rings Closing/M7 (Figure 6) occurs (n49) before Revert To Neck Ring/Blank 
Interference/ml 9 (n51), the collision Neck Rings Collide With Blank Mold/c6 
5 will occur. At nlOO (Figure 5), which is some time (sl3) after the opening of the 
neck rings (M8) at n45/ the invert is displaced back to its original position 
(Revert/M17). Revert is complete at n99/e34/n53. Revert has three submotions 
1. starting at nl00/e33/n48 there is Revert Clears Interference With Blow 
Head/ml 7 which ends at n47/el2/n52, 2. following Revert Clears Interference 

10 With Blow Head, there is Revert To Neck Ring/ Blank Interference/ml 9 which 
ends at n51/el3/n54 when 3. Revert Completion/ m20 operates ending at 
n53/e34/n99. At n50 which is some time (sl4) after nlOO, Neck Rings 
Closing/M7 operates until n49. If the neck rings are not closed before revert 
reaches its initial location of interference with the blank mold (time n49 vs. n51), 

1 5 the collision Neck Rings Collide With Blank Mold7c6 will occur. 

At time nl02 (Figure 5) which is some time (s23) after n23, the movement 
Blow Head On/Ml 8 (Figure 6) takes place finishing at nl01/e36/n59. This is a 
two stage displacement beginning with Blow Head To Interference With 
Revert/m22 which begins at nl02/e35/n58 and ends at n57. In the event that 

20 Revert Clears Interference With Blow Head does not occur before Blow Head To 
Interference With Revert, Revert Collides With Blow Head/c8 will occur (n57 vs. 
n47). The last portion of the blow head displacement is Blow head On 
Completion/m23 which begins at n57/el4/n60 and ends at n59. 

At n56 Blanks Close/M9 (Figure 6) begins and continues until n55. If the 

25 completion of Revert/Ml 7 at n99 does not preceed the start of Blanks Close/M9 
at n56 the Revert Collides With Blank Molds/c7 collision will occur. At nl6 
which is some time (s6) after n99, the Plunger To Loading Position/Ml 
displacement takes place ending at nl 5. 
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n30/el7/n66 (Figure 7) is the beginning of Mold Contact/d5 (Figure 8) 
which ends at n65/el8/n68 and Final Blow/p2 which ends at n63. n30/el l/nl65 
is also the end of Vacuum Blow Lead/dl2 which begins at nl66/e77/nl68. Also 
beginning at nl68 is Vacuum Blow/p5 which ends at nl67 which is some time 
5 (s29) before n68/el8/n65 (Figure 8) which is the end of Mold Contact/d5. Both 
Vacuum Blow Lead/dl2 and Vacuum Blow/p5 begin at nl68/e77/nl66 (Figure 
6) which is some time (s9) after n97. At n91 which is some time (s27) following 
the end of Blow Molds Close/Ml 6 at n97, Blow Mold Cooling/p3 (Figure 8) 
begins continuing until n90 which is some time (s30) before the end 
10 (n65/el8/n68) of Mold Contact/d5. Additionally Blow Mold Precooling/dl 1 
(Figure 6) begins at the same time n91/e74/nl62 and continues until 
nl61/e75/n30/el6/nol which is also the end of Reheat/d4. Finish Cooling/p6 
(Figure 7) begins at nl70 which is some time (s31) after the end of the end of 
Blow Head On/M18 at nlOl and ends at nl69. 
15 At nl04 (Figure 7) which is some time (s32) after the end of Finish 

Cooling/p6 at nl69, Blow Head Up/M19 begins ending at nl03/e38/n73. This 
movement can be broken down into a number of submotions: 1. Blow Head Up 
To End Final Blow/m41 which begins at nl04/e76/nl64 and ends at nl63 which 
is some time (s20) ahead of n63 (the end of Final Blow/p2), 2. Blow Head Clears 
20 Interferencel With Takeout/m25 which begins at nl63/e37/n72 and ends at n71, 
3. Blow Head Up Clears Interference2 With Takeout/m7 which begins at 
n71/e21/n95 and ends at n92, 4. Blow Head Up Clears Interference3 With 
Takeout/m8 which begins at n92/e5/n96 and ends at n94 (Figure 8), and 5. Blow 
Head Up Completion/m26 which begins at n94/e6/n74 and ends at n73. 
25 Tongs Open/MP 12 (Figure 6) is complete at n86 (of the prior cycle t5) and 

some time (s28) thereafter, at nll9, Kickback (takeout ready position)/M22 
begins and ends at nl 18. At nl06, which is some time (s24) after nl 18, Takeout 
In/M20 begins ending at nl05. Takeout movement has a number of sub 
movements: 1. Takeout In To Interferencel With Blow Head/m27 which begins 
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at nl06/e39/n76 and ends at n75, 2. Takeout In To Interference2 With 
Blowhead/m9 which begins at n75/e22/nll7 and ends at nil 6, 3. Takeout To 
Interferences With Blowhead/mlO which begins at nll6/el9/nl32 and ends at 
nl31 and 4. Takeout In Completion/m28 which begins at nl31/e20/n78 and ends 
5 at n77/e40/nl05 (Figure 8). A number of collisions are identified: 1. Blow Head 
Collides With Takeout/c9 which will occur if n75 occurs before n71, 2. Blow 
Head Collides With Takeout/cl4 if nil 6 occurs before n92, and 3. Blow Head 
Collides With Takeout/cl5 (Figure 8)if nl3 1 occurs before n94. At n80 which is 
some time (si 8) after nl05 (the end of Takeout In/M20) Tongs Close/Mil 

10 finishing at n79/e51/nl20. Beginning at n68 and ending at n67/e50/nl22 Blank 
Molds Open/MlO. This motion has a number of submotions: 1. Molds Open To 
Release Point/ m29 which begins at n68/e49/nl21 and ends at nl20/e4/n64, 2. 
Molds Open To Clear Ware/m6 which begins at n64 and ends at nl30/e48/nl23, 
and 3. Molds Open Completion/m3 1 which begins at nl23 and ends at 

15 nl22/e50/n67. At nl08, which is some time (sl9) after n79, the end of Tongs 
Close/Mil, Takeout Out/M21 takes place ending at nl07 (Figure 9). This 
motion also has a number of sub motions: 1. Takeout Out Through 
Interferencel/ml3 which begins at nI08/e41/n!38 and ends at nl33, 2. Takeout 
Clears Ware From Mold/m30 (Figure 9) which begins at nl33/e57/n82 and ends 

20 at n8 1 , 3 . Takeout Out Through Interference2/m24 which begins at n8 1 /e23/n 1 3 5 
and ends at nl3, 4. Takeout Out Through Interference3/m36 which begins at 
nl3/e58/nl37 and ands at nl36, and 5. Takeout Out Completion/m37 which starts 
at nl36/e59/n88 and ends at n87/e42/nl07. The collision Molds Collide With 
Takeout/cl 1 will occur if n82 occurs before nl30. 

25 Finally at the end of Takeout Out/M21 (nl07/e46/nll5)Hanging Dead 

Plate Cooling/d6 takes place until nl 14. At n 174 which is some time (si 2) after 
nl07, Dead Plate High/p8 takes place lasting until nl8. Some time (s33) 
thereafter, at n84/e78/nl76/e47/nl 14 Total Process/dl3 ends, Hanging Dead 
Plate Cooling/d6 ends and Tongs Open/Ml 2 ending at n83. 
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While, for illustrative purposes, one specific blow and blow machine 
configuration has been described, it should be understood that there are a number 
of operating configurations which machine users use including blow and blow 
and press and blow and for each users have developed many unique processes 

5 that would vary slightly one from another. A person skilled in this art, with an 
understanding of the illustrated configuration, should be able to define a 
constraint diagram for his actual configuration. 

The next step is to convert this network constraint diagram into a 
representation that is ideal for automated formulation and solution of the schedule 

1 0 synthesis and analysis problems by a computer. A matrix algebraic representation 
of the network constraint model is utilized in the preferred embodiment but other 
forms of mathematical representations can be used. The Branch Incidence 
Matrix, F may be formed as follows: 

1. Number the branches in the network constraint diagram (NCD) from 1 to 
1 5 M b where, M b is the total number of network branches. The ordering of the 

assigned branch numbers is arbitrary. 

2. Number the nodes in the NCD from 1 to N n where N n is the total number 
of network nodes. The ordering of the assigned node numbers is arbitrary. 

3. Form the first row of an M b row by N n column matrix F by entering a value 
20 of 1 (positive one) in the column corresponding to the source node for the 

first branch, a value of -1 (negative one) in the column corresponding to the 
destination node for the first branch, and zeroes in all the other columns. 

4. Create the second through M b row of F by repeating the procedure described 
in Step 3 for the second, third on up to the M b branch in the network. 

25 

The result will be a matrix, F with M b row by N n columns which is almost 
entirely filled with zeros, except for one entry of 1 and one entry of -1 in each 
row. 
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To provide a concrete example, the NCD for a simple network model is 
shown in Figure 10. The network has M b =7 branches, and N„ =6 nodes. The 
Branch Incidence Matrix, F for this network will thus have 7 rows and 6 columns. 
For this model, utilizing the branch and node numbers indicated in Figure 3, F 
5 will then be given by: 



f 1 — 1 0 0 0 0] 

| 0 1 0 0 -1 0 | 

| 1 0 -1 0 0 o I 

10 F= I 0 0 0 0 1 -1 i Equation 1 

i o o i—i o o ! 

| 0 0 0 1 0 -1 I 

L o i o-i o oj 

1 5 Each branch, i, in the network constraint model represents a pair of 

relationship of the form: 

t destination,! " t sourcej — § m ax,i Equation 2 

t destination^ - t source,! ^ S m in,i Equation 3 



20 



Where: 



t destination,; = time assigned to the destination node of the i th branch 
t source i = time assigned to the source node of the i th branch 



25 5 max ,i = maximum allowable branch duration for the i th branch 
S min ,i = minimum allowable branch duration for thei th branch 
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Define the vector t of nodal times, where the j n element oft is the time 
assigned to the j* network node. Denoting the i th row of the branch incidence 
matrix F by F„ Equation 2 and Equation 3 can be rewritten as follows: 

5 -Fit < 8 max ,i Equation 4 

-Fjt > 8 min>1 Equation 5 

This results from the fact that that the matrix multiplication of the i th row 
10 of the constraint matrix,Fj, with the nodal time vector, t, selects only the source 
and destination node times, because all other entries in the row are zero. In 
accordance with conventional practice a value of positive one is assigned to the 
element corresponding to a source node and a value of negative one is assigned to 
the destination node. 

1 5 Since Equation 3 and Equation 4 hold for each branch in the network the 

Fundamental Matrix Constraint Equations can be written as follows: 

-Ft < 5 max Equation 6 
-Ft > 5 min Equation 7 

20 

For branches that have no upper bound on their duration, S,„„ v ,/ is set to 
positive infinity. Similarly for branches which have no lower bound set on their 
duration, 8 min j is set to negative infinity. For branches which must exactly match a 
target value, both the upper and lower bound are set equal to the target value, 

25 Starget,i- 
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The fundamental matrix constraint equations (Equation 6 and Equation 7) 
are augmented to account for three types of additional constraints. These three 
additional types of constraints are: 

5 1 . The branch duration for all cycle branches must be equal. This is required to 
enforce a uniform cycle period on the whole system. 

2. The branch duration for each submotion branch must remain a constant 
fraction of the branch duration for its corresponding main motion branch. 

10 

3. The absolute event time for one node in the network must be set to a desired 
reference value (typically zero). 

These requirements can be expressed in terms of the previously defined 
1 5 branch incidence matrix, F as follows: 

Each cycle branch must have a duration equal to the cycle period, T, and 
therefore, regardless of the particular value of the cycle period all N t cycle 
branches must have the same branch duration. Denote the branch numbers 
20 corresponding to the cycle branches by the set {ii,i 2 ,-iN t } • The duration of the k th 
cycle branch can then be expressed as: 

-Fj k t = 5j k Equation 8 
25 Where Fi, represents the k row of the Branch Incidence Matrix, F. 
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The uniform branch duration constraint is then enforced by setting each of 
the cycle branch durations to be equal to the duration of the first cycle branch as 
follows: 

r -F i2 i r -Fij i 

1 -Fi 3 l t = l -Fj, l t Equation 9 

L-fJ L-fJ 

Which can be rearranged to obtain: 

r -F i2+ F tl i r o i 

I ^3 + ^, I t = I 0 I Equation 10 
L -Fj N + Fs, J L 0 J 

Defining the matrix A t to represent the left hand side of Equation 10 so that: 

r -F i2+ F ;i i 

A t = I "^'3 + F 'i I Equation 1 1 

i i i 

L -F iN + Fi, J 

Equation 10 can be rewritten more compactly as: 

A t t = 0 Equation 12 

If the branch duration for a main motion branch is varied, then the 
submotion branches (if any) associated with this branch should be 
proportionally rescaled. 
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20 



25 



To represent this set of auxiliary constraints some notation must first be 
defined. Denote the branch numbers corresponding to the main motion 
branches (including only branches that have associated submotion branches) by 
the set {M,,M 2 ,-M Nm } where Nm is the total number of main motion branches 
5 that have associated submotion branches. Denote the submotion branches 

associated with the k th main motion branch by the set {m kl ,m k2 ,-m kNk }, where 
N k is the total number of submotion branches associated with the k th main 
motion branch. Each submotion branch duration represents a fixed fraction of 
the associated main branch duration. Let a kj denote this fixed fraction for the j th 
1 0 submotion branch associated with the k th main motion branch. 

The required set of constraints associated with the kth main motion 
branch can now be represented by the equation: 

r F mkl -a kl F Mk l r 0 ] 
15 I F mk2 -a k2 F Mk l t== | 0 I Equation 13 



L F m - a kN F Mk J L 0 J 
Defining the matrix A Mk to represent the left hand side of Equation 13 so that: 

F F mk] - a k iF Mk 1 

A Mk = ' Fm k2~ ak2pM k I Equation 14 
L Fm kN k ~ akNFM J 
we can rewrite Equation 13 more compactly as: 

A M| t = 0 Equation 15 
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Further, defining the matrix A m as: 



r 

I 



the complete set 
equation: 

10 

A m t = 0 Equation 17 

One reference node in the network is selected, and the absolute time at 
which this event is to occur is set to zero. Denoting the node number of the 
1 5 reference node as k, this constraint can be expressed as: 

A z t = 0 Equation 18 

where the kth element of the row vector A z has a value of 1 and all other 
20 elements are zero. 



A M , I 
Am 2 I 

A M3 I Equation 16 

A M NM -J 

of auxiliary submotion constraints can be expressed by the 



Finally, the Augmented Constraint Matrix A is defined by: 

r -f, i 

25 A= I A t I Equation 19 

I A m I 

L A z J 
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where F,- is the reduced branch incidence matrix formed by eliminating all of 
the now redundant rows in F. Specifically, the rows corresponding to all of 
the submotion branches, and all but the first cyclic branch are removed from F 
5 to form ¥ r . The length N*+N,+Na/ +1 vectors b mi „ and b nmx are defined by: 



b min = I I Equation 20 

L o J 

10 

T §max 1 

b max = I I Equation 21 

L o J 



1 5 The complete set of network constraints can then be expressed by - 

combining Equation 6, Equation 7, Equation 12, Equation 17, Equation 20 and 
Equation 21 into the single set of augmented constraint equations: 

At < b max Equation 22 

20 

At > b min Equation 23 



In the right hand side of Equation 22 and Equation 23 the final 
N t +N M +1 elements of the vectors b max and b min , respectively, are always 
25 identically zero. The objective is to find a set of event times (schedule) which 
satisfy all the required network constraints. In general, there will be more than 
one, in fact, infinitely many, schedules that will fit the network constraints. A 
constrained optimization based methodology is therefore employed to select 
the most desirable schedule from the many available possibilities. The general 
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approach can be tailored to fit a variety of practical problems by making 
judicious choices of the optimization criterion. Some informally stated 
examples of optimization criteria that are of practical interest include: 

5 1 . Minimizing the cycle period with specified thermal process durations 

2. Maximizing particular thermal process durations, for example reheat, within 
a fixed cycle period. 

3. Minimizing wear and tear by slowing down mechanisms as much as possible 
with a fixed cycle period and specified set of thermal process durations. 

Optimal schedules utilizing such criteria are readily obtained using the 
new methodology which has been developed. 

In terms of the matrix algebraic model representation, described 
previously, the general problem to be solved is to find a length N n vector t of 
nodal times which satisfies: 

minimize f(t) Equation 24 
20 subject to the constraints: 



10 



15 



The scalar function f, referred to as the objective function, specifies the 
criterion for distinguishing the most desirable of the many possible solutions to 
the problem. This is known as a Constrained Optimization Problem (as 
opposed to an Unconstrained Optimization Problem) because we are seeking an 
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optimal solution but are limiting the set of possible solutions to those that 
satisfy a specified set of constraints. In this case, the constraints are expressed 
as a set of linear inequalities. 

A wide variety of practical criteria can be expressed in terms of a 
5 quadratic objective function of the form (actually, the constant term,f 0 , is not 
strictly required since it has no effect on the location of the system's minima 
and maxima. It is only retained here because it later allows the value of the 
objective function to be given a more obvious interpretation as the distance of 
the actual branch durations from the desired target values. 

10 

f(t) = Ya^Ht+Ct+fo Equation 25 

As will be detailed subsequently, the essential machine scheduling 
problems can in fact be expressed using a quadratic objective function of the 
1 5 form given in Equation 25. 

An optimization problem that has this combination of a quadratic 
objective function, and linear constraints is known as a Quadratic Programming 
Problem. A large variety of fast and reliable numerical algorithms exist for 
solving Quadratic Programming Problems. In some practical cases (for 
20 example minimizing the cycle period) the optimization criteria can be 
expressed using a linear objective function of the form of given: 

f(t) = Ct + f 0 Equation 26 

25 This combination of a linear objective function with linear constraints is 

known as a Linear Programming Problem. Linear Programming problems, in 
many cases, can be solved with less computational effort and therefore even 
faster than Quadratic Programming Problems, but the Quadratic Programming 
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Solver used to save linear and quadratic objective functions is the most 

economical solution. 

The basic idea of the GTSSM (General Target Schedule Synthesis 

Methodology) is to assign a target value for the duration of each branch in the 
5 network. These target values represent the ideal set of values the user would 

like to attain for all of the branch durations. Because of the many network 

constraints that must also be satisfied, it may not in fact be possible to obtain 

all of the target branch duration values. The GTSSM therefore finds a schedule 

that matches the target values as closely as possible. 
10 The GTSSM achieves it ability to provide a single approach to a variety 

of problems through the use of four major features: 

1. Quadratic Objective Function - A quadratic objective function makes 
mathematically precise the notion of a schedule being as close to the target 

1 5 value as possible. 

2. Hard Limits - Hard upper and lower limits can be imposed on the allowable 
durations for each network branch. 

20 3. Locking - Durations of specified branches can be locked so that they are 
achieved exactly in the resulting schedule. 

4. QP Solver - Use of a robust QP (Quadratic Programming) numerical solver. 
Each of the above features will now be described in further detail. 
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The intuitive notion of a schedule being close to the target value must be 
made mathematically precise in order to implement an automated numerical 
solution. For this purpose define the objective function, f(t), as follows: 



5 f(t) = Z,l;( w .^( t )-^,)) 2 Equation 27 

where: 

= a constant which weights the importance of the deviation between the 
target and actual duration for the i th network branch 

= duration of the i th network branch as a function of t, the length N„ 
vector of nodal event times (schedule) 

5 t = target duration for the i th network branch 
15 

N b = total number of network branches 

Thus, the distance, from the target is expressed as the weighted, sum of 
the squared deviations between the target and actual branch durations. It is 
20 noted that for the two or three-dimensional case (N b =2 or N b =3) and w ; =l , 
Equation 27 expresses the familiar Euclidean Distance Formula. 

Noting that the branch duration for the i th branch duration can be 
expressed in terms of the i th row of the branch incidence matrix as: 

25 5j = -Fit, 



10 



5i(t) 
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Equation 27 can be expressed in terms of the previously defined Matrix 
Algebraic System Model as: 

f(t) = (W(Ft+5)) T (W(Ft+5)) Equation 28 

5 

where: 

W= weighting matrix 
5= vector of target branch durations 
1 0 F= branch incidence matrix 

T= length N n vector of nodal event times (schedule) 
superscript T = matrix transpose 

Following routine algebraic manipulation Equation 28 can be rewritten as: 

15 

f(t) = t T F T W T WFt + 28 t W T WFt + 5 t T W T W5 t Equation 29 

Equation 29 can then be expressed in the standard form given in Equation 25 
for a quadratic objective function: 

20 

f(t) = ^Ht+Ct+fo Equation 30 

where: 

25 H= 2F T W T WF 

C= 25 t W T WF 
f 0 = §t T W T W5 t 
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Some flexibility is available in the definition of the elements of the 
diagonal weighting matrix. The simplest alternative is to set each of the branch 
weights Wj to a value of 1 (one) so that W becomes the identity matrix. This 
gives equal weight to the absolute error (deviation) between the desired and 
5 target duration values for all of the network branches. Although in some cases 
the absolute error approach may be appropriate, it is more commonly the case 
that we are concerned with the relative error, in which the error for each branch 
is normalized by its typical duration. Thus, with the relative error approach a 
deviation of 1 millisecond for a branch whose typical duration is 10 
1 0 milliseconds is considered to be of the same significance as a deviation of 1 

second for a branch whose typical duration is 10 seconds. For the relative error 
approach we thus define the weighting matrix W by: 

15 T l/(5 higV 5 l0Wl ) 0 0 

W= ' ® l/(8 hjgh2 -§i OW2 ) 0 

i 0 0 

L o oo 

20 where: 

Shighi = The high scale value for the i th network branch 
Siowi = The low scale value for the i th network branch 

25 It is often useful to have the ability to limit the allowable ranges of 

particular branch durations. Example situations for requiring this capability 
might include mechanisms that have a lower bound on their motion duration, 
and process steps that have lower and/or upper bounds on their duration. These 



0 1 

0 I Equation 3 1 

0 I 

l/(Shigh N -Slow N ) J 
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bounds are set in the GTSSM by assigning appropriate values to the elements 
in the b min and b max vectors forming the right hand sides of the matrix constraint 
relationships given by Equation 24. 

In some cases it is desirable to specify that particular branch durations 
5 are to be exactly equal to the target values. This will be referred to as locking 
the target value. For example, in some cases it is necessary to lock the duration 
of the cycle branches because the cycle period of upstream equipment such as 
the feeder can not be readily adjusted. This capability is implemented in the 
GTSSM by setting the value of the appropriate elements of the upper and lower 

1 0 bounds (vectors b min and b max forming the right hand sides of the matrix 

constraint relationships given by Equation 24) both equal to the target value. 
The matrix H should be positive definite). To avoid complications with this 
numerical problem, small weights can be assigned even to branches whose 
duration's are not of interest, or a solver can be used that specifically treats the 

1 5 case where H is only positive semidefinite. 

Based upon either prior experience or specific tests, the desired 
durations of all of thermal forming process steps (reheat, final blow, etc.) may 
be known and the bottle maker may not wish these values to vary. With the 
cycle period branch duration unlocked, all of thermal-forming related branch 

20 durations locked, and mechanism motion branches duration locked at the value 
corresponding to the fastest possible mechanism travel durations, the target 
duration for the cycle period could be set to zero (indicating that it is desired to 
have it as short as possible). The QP solver will then find the schedule with the 
shortest possible cycle period consistent with all of the network constraints 

25 (These constraints include the locked thermal process durations and mechanism 
motion durations along with the requirement for collision avoidance, proper 
sequencing, etc.). 

It is quite possible that a particular schedule may achieve the required 
cycle period and desired set of thermal formal process durations but that it 
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requires moving some mechanisms faster than is strictly necessary to achieve 
these objectives. It may be desirable instead to operate the mechanisms only as 
fast as is absolutely necessary to achieve the other desired objectives. This 
would reduce the average and peak current to the servomotors (and associated 
5 motor heating) and perhaps otherwise reduce general wear and tear on the 
system. To do this, the cycle period and other thermal process branch 
durations would be locked at their desired values. All of the motion branch 
durations would be unlocked, and their target values set to a relatively large 
value. The QP solver would then have the freedom to speed up the mechanism 
10 if required to meet the constraints on cycle period and thermal process 

durations, but otherwise it would increase the motion durations as much as 
possible. 

When the desired target values can not be exactly achieved, the user can 
be provided with some indication of which bounds must be relaxed in order to 

15 more closely achieve the desired objective. This can be done by examining the 
Lagrange Multiplier Values computed at the location of the optimum schedule. 
The Lagrange Multipliers can be given the interpretation of being the partial 
derivatives of the objective function with respect to the elements in the b min and 
bmax vectors forming the right hand sides of the matrix constraint relationships 

20 given by Equation 24. Thus non-zero values for a particular Lagrange 

Multiplier indicates that the objective function would be either increased or 
decreased (depending on the algebraic sign of the Lagrange Multiplier) by 
changing the value of the associated element of the b min and b max vectors. Such 
constraints are said to be active. Other constraints whose Lagrange Multiplier 

25 values are zero are said to be inactive. By appropriately displaying to the user 
the active constraints ranked by the relative magnitude of their Lagrange 
Multiplier values the user would be informed as to which bounds are imposing 
the greatest limitation on achieving the desired results. Further, the sign of the 
Lagrange multiplier could be used to determine, and subsequently display to 
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the user, whether the target value (in the case of a locked branch) should be 
increased or decreased to further improve the ability to achieve the target 
values of the unlocked branches. Most constrained optimization algorithms 
provide for the capability to compute Lagrange Multiplier values (or already 
5 compute them as part of their normal operation), so this additional information 
could be utilized to provide further guidance to the user if desired. 

If the user overly constrains the system there may be no feasible solution 
to the QP problem that has been posed. In such a case it is important to 
recognize that the problem is infeasible, and to relax the bounds enough to 

1 0 allow a feasible solution to occur. QP solvers typically are able to recognize 

that there is no feasible solution and return an appropriate flag. This flag can be 
used by the software that implements the GTSSM to trigger a prompt to the 
user to relax any constraints as much as possible. 

The MAR (Matrix Algebraic Representation) also allows a proposed 

1 5 schedule to be analyzed in order to discover any potentially damaging or 
undesirable constraint violations. This capability provides a mechanism for 
performing intelligent input qualification on user requested changes to event 
times which goes well beyond conventional high and low bound checking. 

The basic purpose of the schedule analysis methodology is to provide 

20 the capability to check a proposed schedule for constraint violations and then to 
report any violations that might be found. The method also allows violations to 
be reported in a way that allows a user to understand the consequences of the 
violation and, to the extent possible, indicates a remedy. 

The actual checking of the constraint violations is computationally quite 

25 simple involving only a single matrix multiplication and subtraction. To fully 
obtain the desired functionality there are some additional considerations that 
must also be taken into account. The additional complexity arises primarily 
from the fact that the can only schedules nodal (event) times for a subset of the 
nodes in the, overall system model. This subset of nodes is referred to as the set 
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of independent nodes. The nodal times for the remaining, dependent nodes, is 
then automatically calculated from the independent nodal times and known, 
fixed branch durations. 

The overall methodology then consists of the following components: 
5 1 .Solving for Dependent Nodal Times 

2 .Detecting Constraint Violations 

3 .Diagnosing and Categorizing Violations 

1 0 Dependent Nodal Times can be solved for in terms of the previosly 

defined constraint set by utilizing the following procedure. 

1 . Form the subset of equality constraints as: 

Aeg t = b eq Equation 32 

15 

By retaining only those rows of A and b (as defined in Equations 19 
and 20 respectively) for which the upper and lower bounds are equal. Note that 
the upper and lower bounds for branches with known, fixed, durations will both 
be set to this known fixed value. The upper and lower bounds of these fixed 

20 duration branches will therefore be equal and the rows of A corresponding to 
these branches along with the auxilliary constraints will thus be retained in A e . 
Typically the branches with known fixed values will be the Motion, Cyclic and 
Simultaneity branches. In order to have a well posed problem, the row 
dimension of A eq must be greater than or equal to the number of dependent 

25 nodal times. It is required that a sufficient number of branches be assigned 
fixed values such that this condition is met. 

2. By reordering the columns of A eq , form the partitioned incidence matrix A p , 
in which the first Njcolumns of A p correspond to the independent nodal 



32 



times. Form the partitioned nodal time vector t p by sorting the columns oft 
to correspond to the new column order in F p , Equation 32 can then be 
rewritten as: 



r t P1 i 

[A P A P ] I — I =b eq Equation 33 



10 3. Rearrange Equation 33 to form the set of linear equations 
A Pd tp D - (beq-A^tp,) = 0 Equation 34 



15 5. Assign current values to the independent nodal event times and elements of 
b eq corresponding to fixed branch duration and solve the overdetermined 
system of Equations 34, for t p . This can be done using standard numerical 

methods available for solving overdetermined sytems of linear equations, e.g. a 
Linear Least Squares Solver. For a consistent set of fixed branch durations and 
20 properly constructed Network Constraint Diagram an exact solution to this 

overdetermined problem may be obtained. That is, a vector of dependent nodal 
times tp D can be found that satisfies Equation 34 without any error. If an exact 

solution can not be found then the user should be notified accordingly so that 
the situation can be remedied. It is noted that the zero reference node should be 
25 included with the independent event times and be consistent with the definition 
provided in Eqation 18. 
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6. The elements of of t p and t p are resorted into their original order 
corresponding to the rows of Equations 22 and 23 to form a vector of times, 
^proposed consituting the proposed schedule. 

5 Once the dependent times have been calculated and a proposed schedule 

is available, actually detecting constraint violations is relatively 
straightforward. Let the proposed schedule be given by the vector of nodal 
times, t pr oposed- From Equation 22 and Equation 23 the conditions to be checked 
are then given by the set of inequalities: 

10 

Atp r0 p 0S ed - b max < 0 Equation 35 
At proposed - b min > 0 Equation 36 

15 if the inequalities given in either Equation 35 or Equation 36 are not 

completely satisfied then the proposed schedule violates at least one constraint. 

Each row in Equation 35 and Equation 36 represents a particular system 
constraint. Explanatory text and a severity level can accordingly be assigned to 

20 each row in these equations. A proposed schedule would then be tested by 

evaluating Equation 35 or Equation 36. The row numbers of any rows that did 
not satisfy the required inequality would then provide an index for recalling 
and displaying corresponding error message text. The severity level could be 
used to sort multiple constraint violations in order of severity, and could also 

25 key for an appropriate color code, or other attribute (flashing) on the graphical 
user interface. 

This assignment of text and severity level can be done in an automatic 
manner. To understand how such an automatic assignment can be performed, 
recall that the rows in Equation 35 or Equation 36 are derived from network 
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branches. The violation implied by each branch type can therefore be an 
attribute that is assigned to the particular branch type and then further specified 
for the particular branch. For example for a collision branch we could 
automatically define the violation text to read "collision occurs between invert 
5 and baffle" this event could also be assigned a severity level for example a 
number between 1 and 10 with 10 being the most severe. The corresponding 
row in Equation 35 or Equation 36 would then inherit these descriptions from 
the branches they descend from. Alternatively, once a Network Constraint 
Diagram was fully defined for a particular forming process, individual 

10 messages could be manually entered or the automatically generated default set 
could be edited and the resulting data stored in a table for each of the finite 
number of constraint violations which might occur. While this manual 
approach could perhaps allow some enhancement to the readability of the 
messages, it could also be error prone, and would have to be updated if any 

15 changes were made to the Network Constraint Diagram. The automatic 
approach is thus preferred. 

In a state of the art control, the operation of one of these 
mechanisms/processes is controlled by turning the mechanism, etc., "on" and 
"off at selected angles within a 360° cycle. The turning "on" of a mechanism 

20 is an event and the turning "off of a mechanism is an event. Figure 12 

illustrates a conventional list of timed events with their angular "on" and "off 
times for an I.S. machine. This list is available from the machine control. 

The unwrapped schedule can be converted to a corresponding wrapped 
schedule using the known cycle period and calculating event angles modulo 

25 360 degrees (event angle = modulo 3 60(un wrapped event time/cycle period) x 
360. To go from a wrapped schedule to an unwrapped schedule the original 
network constraint diagram is augmented with a new set of directed branches 
called unwrapping branches. The subgraph formed from the unwrapping 
branches along with the motion and sequence branches and any nodes that are 
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incident on these branches will be referred to as the cycle unwrapping graph. 
An example of a Cycle Unwrapping Graph is illustrated in Figure 12 which 
shows a press and blow cycle. The CUG (Cycle Unwrapping Graph) is created 
so as to have the following properties. 

5 

Property 1 . The CUG is a connected graph 

Property 2. The nodes of the CUG are exactly the set of all the source and 
destination nodes for all of the motion and process branches in the NCD 
10 (Network Constraint Diagram). This means that every "on" and "off angle of 
the timing drum (sequencer) is represented on the graph. 

Property 3. Every branch in the CUG is part of a cycle (path from an event to 
the next periodic repetition of that event). For example, the lowermost line of 

1 5 figure 1 1 proceeds : M 1 20 (Tongs Open), MHO (Tongs Close), M2 1 0 

(Takeout Out) and Ml 20. Similarly, the next line up proceeds M210 (Takeout 
Out),M220 (Kickback), M200 (Takeout In), and M210. The next line up 
proceeds, M190 (Blow Head Up), M180 Blow Head Down), p2 (Final Blow), 
and Ml 90. The next line up proceeds MP 1000 (Blow Molds Open), M240 

20 (Blow Molds Closed), M1000. The next line proceeds MP 100 (Plunger To 
Loading Position), M230 (Pressing), M80 (Neck Rings Open), M70 (Neck 
Rings Closing) and M100. The next line proceeds MP90 (Blanks Close), 
M230, M40 (Plunger To Invert Position), M60 (Invert), M70 (Revert), M90 
and the first line from the left proceeds MP 150 (Baffle Up), Ml 40 (Baffle 

25 Down), M230.M150 (Baffle Up). 



Property 4. The branches incident on the nodes of the CUG either fan in or fan 
out but not both. That is, if there is more than one branch directed towards a 
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given node then there is only one branch leaving that node (fan in). If there is 
more than one branch leaving a given node then there is exactly one node 
entering it (fan out) 

The above Properties imply that the CUG also has the following 
5 additional properties: 

Property 5. Any sequence of three nodes that can be traversed by two 
interconnecting branches following in the branch direction, will be part of at 
least one common cycle. So each node in the CUG is between two other events 
10 in a cyclic sequence. 

Property 6. Since each branch in the CUG is part of a cycle it must be less than 
one period long. 

The problem is then solved in a series of steps , which include checks 
15 that the input data is properly ordered to provide a valid solution. 

1 . Form the branch incidence matrix for the CUG. 

2. Divide the nodes of the CUG into two sets: The independent nodes whose 
20 values are given in the input set of wrapped event angles, and the remaining 

dependent nodes whose event angles are as yet unknown. For a well posed 
problem all of the dependent nodes must be connected to an independent 
node by a branch whose duration is known. 

25 3. Assign the known input event angles to the independent nodes in the CUG 
to which they correspond. 
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4. Determine the event angles for the dependent event angles using 



(h). = mod (>-<). ± 




di-i 




cycle 



*360 ,360 



J 



Equation 37 



Where: 



5 0 is the event angle to be calculated for ith dependent node 

(h) . is the dependent node connected to node i through a branch 



with known temporal duration . 

The algebraic sign in Equation 37 is chosen as positive when the dependent 
node is downstream from the independent node and negative otherwise. 
10 5. Assign the event angle for any periodic repetition nodes to equal the value 

of the node which it replicates (the node to which it is connected to by a 

cyclic branch in the NCD). 

6. Check that all event angles are in the correct cyclic order. This is done by 
checking that the event angle assigned to each node is between the value of 

15 any pairing of its immediate upstream and downstream adjacent nodes. 

7. Find the angular branch durations for all of the branches in the CUG using: 

5 = mod(- F@,360) Equation 3 8 

20 Where: 

F is the branch incidence matrix for the CUG 
© is the vector of nodal event angles in the CUG 
8 is the vector of angular branch durations in the CUG 
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8. Convert 5 , the vector of angular branch durations to a vector d of 
temporal duration using: 

d = 5/360* T Equation 39 

5 

where T is the cycle period. 

9. Solve for the unwrapped nodal event times by solving, using standard 
numerical methods, the possibly overdetermined system: 

10 - F r t = 6 Equation 40 

Where F r is the branch incidence matrix of the CUG with the column 
corresponding to the zero reference node deleted. (The choice of zero reference 
node is arbitrary but should be consistent with that of the NCD.) Although the 
1 5 above system is overdetermined the least squares solution will in fact have zero 
error because 6 is in the column space of F r . This should be verified to identify 
any computational problems. 

10. The independent nodes in the NCD are assigned values using the 

20 corresponding unwrapped event times that are determined from equation 4. 
The dependent nodes in the NCD can then be determined as previously 
described. 

Figure 13 is a block diagram illustrating the making of the analytical tool 
25 (Tool). The first thing is to Define A Network Constraint Diagram For A Bottle 
Forming Process In An I.S. Machine 60 (an unwrapped cycle following the 
formation of the gob, its delivery to the blank station, the transfer of a parison 
from the blank station to the blow station and the removal of a formed bottle from 
the blow station). Then Translate The Network Constraint Diagram Into A Data 
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Table 61 which is a compilation of the essential data in the constraint diagram 
and includes lists of all branches and for each, identifies its name, type, start node 
number, end node number, etc.. This can be done manually by examining a 
drawing of the constraint diagram or a computer can do this from a computerized 
5 drawing of the constraint diagram. The next step is to Translate The Data Table 
Into A Mathematical Representation 62 which is then transformed into a 
Computerized Model 64 for evaluation. In the illustrated embodiment the 
computerized model is based on matrix mathematics but other mathematical 
approaches could be used. Depending on the nature of the bottle and the process 

10 used to make the bottle (press and blow, blow and blow, for example), different 
data tables may be required. As shown, a number of data tables (Data Table "N") 
may be stored in memory and input as desired. 

Figure 14 illustrates the portion of the Computerized Model 64 which will 
Unwrap 360° Machine Cycle Event Angles Into Bottle Forming Process Event 

15 Times 66 (Unwrapping). It receives, as inputs, either as downloaded data from a 
machine control 42 (fig. 1) or the like, or via manual entry at a console, terminal 
or the like, the Event Angles, the Machine Cycle Time (the cycle time for a 360° 
timing drum for an I.S. machine) and the Motion Durations (the "M" (capital M) 
durations for the displaceable mechanisms) and outputs Event Times in the bottle 

20 forming process. It should be noted that while the event angles and machine 
cycle time are normally available data from an existing job file, motion durations 
would have to be defined for the job. 

Figure 15 illustrates the use of this Computerized Model 64 to Analyze An 
Unwrapped Schedule Re Constraints 68 (Constraint Violation Avoidance). With 

25 Event Times, Machine Cycle Time, Motion Durations, Submotion Durations (the 
"m" (small m) motion durations for the displaceable mechanisms), the Collision 
Branch Lower Limits, the Sequence Branch Lower Limits, and the Thermal 
Forming Process "N" Limits as inputs, the computerized model can determine 
whether There Will Be A Constraint Violation? 70. While the word input means 
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that the inputted information is made available, it might be available form 
different sources. For example, the event angles and machine cycle time could be 
available from an existing job file, whereas the rest of the inputs would be 
inputted at the time when the data table is inputted into Translate The Data Table 
5 Into A Mathematical Representation 62. 

Whenever an input can have a range of values that may be selected by the 
operator, such input will include the upper and lower limits of that input, and a 
choice as to whether the setting is to be locked at a specific value or unlocked to 
permit its location somewhere within the limits. Nominally the lower limits for 

10 the Collision and Sequence Branches can be set at zero or at a selected margin of 
error and this can be locked out from the operator or the operator may be given 
access to these inputs so that the operator can define any desired lower limits. 
One constraint violation would be a schedule that would result in something 
happening in the wrong sequence. Another would be a schedule that would result 

15 in a collision. Either of these constraint violations could be determined without 
Thermal Forming Process Duration "N" Limits. With this additional input(s), the 
unwrapped schedule could be evaluated to determine whether one or more of the 
Thermal Forming Process Durations will be either too short or too long and 
thereby violate one or more of the thermal forming process constraints. These 

20 inputs and ouputs, as well as inputs and outputs in latter discussed embodiments, 
could be available for viewing on any suitable screen. 

If either inquiry obtains an affirmative answer, the control will Operate 
Alarm and/or Reject Inputs 74 and Output Constraint Violation(s) 76. If neither 
inquiry is answered in the affirmative the control can Output The Calculated 

25 Margins 78 to give the operator some idea as to how tight the schedule is and then 
Wrap The Event Times Into Event Angles And Print The Event Angles and New 
Machine Cycle Time 79. "Print" is intended to mean the presentation of data in 
either operator readable form as an output presented on a screen or a document or 
machine readable form so that the machine control can automatically operate on 
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the data such as by resetting the machine with the new event angle and machine 
cycle time. 

In one mode, an I.S. machine may be running and the operator may want 
to change one or more of the event angles in the 360° timing drum. A particular 
5 job is being run and basic data for that job (the durations and limits) has already 
been inputted into the control. This data along with the machine cycle time can 
be downloaded from the machine control. The Event Angles including any 
proposed Event Angle change can be downloaded to the unwrapper 66 so that 
Event Times can be defined. In another mode, an operator may have a record 

10 (Event Angles and Machine Cycle Time) of a job that was run previously and 
want to evaluate some changes before he starts the job. 

In a conventional I.S. machine which has a number of mechanisms that are 
operated via pneumatic cylinders, Motion Durations and Submotion Durations 
may have to be empirically defined, as with high-speed cameras. Where 

15 interferences involve actuators which are displaced pursuant to motion profiles, 
the submotion zones can either be empirically defined or they may be 
mathematically determined. 

Figure 16 illustrates the use of this computerized model to monitor the 
thermal forming process durations (Thermal Forming Process Durations). With 

20 Event Times, Motion Durations, Submotion Durations and Machine Cycle Time 
known or as inputs, the Computerized Model 64 will Analyze An Unwrapped 
Schedule Re Thermal Forming Process Durations 80. and then the Computerized 
Model 64 will Output The Thermal Forming Process Durations 82. The operator 
can, accordingly, at any time, see the Thermal Forming Process Durations and 

25 based on his experience, make changes to the 360° Event Angles and Machine 
Cycle Time. With the additional input of Thermal Forming Process Duration "N" 
Limits, the computer model can also Output Thermal Forming Process Duration 
"N" Margins 81 so that the operator can see where the time of any process is 
relative to its allowable time window. 
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Figure 17 illustrates the use of the computerized model to define, for an 
existing machine set up, the optimized cycle time (Optimized Cycle Time) and 
the optimized Event Angles for that schedule. With Motion Durations, 
Submotion Durations, Collision Branch Lower Limits, Sequence Branch Lower 
5 Limits, Event Times, Machine Cycle Time, and Optimized Machine Cycle 
Time/Target/Lock Status known or as inputs to the Optimize Unwrapped 
Schedule For Minimum Cycle Time 82, the Computerized Model 64 will 
determine whether There Is A Feasible Schedule? 83. If not the model will 
Reject The Inputs 85. The Machine Cycle Time and the Event Times may be 

10 supplied from the unwrapper and the Optimized Machine Cycle Time may be 
inputted by the operator. The Event Times and Machine Cycle Time are only 
required to determine the thermal forming durations so that these valves could be 
locked before doing the optimization. Equivalent inputs would be the Thermal 
Forming Durations. The operator can set the Optimized Machine Cycle Time 

15 Target to zero with an unlocked status and the Computerized Model will try to 
optimize the proposed schedule at the lowest possible cycle time. In the event the 
operator decides that rather than reduce the machine cycle time from the current 
Machine Cycle Time to the fastest Machine Cycle Time, he would prefer to 
reduce the cycle time to some machine cycle time therebetween. He can set the 

20 Optimized Machine Cycle Time Target at a time intermediate the Machine Cycle 
Time and the fastest machine cycle time with a locked status. If there is a feasible 
schedule, the model will Wrap Optimized Event Times Into Event Angles 84 and 
Print The Event Angles And The New Machine Cycle Time 86 for the schedule 
cycle so that it will be available for input into the machine controller portion of 

25 the control. 

Figure 18 illustrates the use of the Computerized Model 64 to tune an 
operating LS. machine in response to operator inputs defining one or more of the 
Thermal Forming Process Durations (Thermal Forming Process Duration "N", 
and associated Target, Limits and Lock Status). With Machine Cycle Time and 



43 



Event Times (or Thermal Forming Process Durations), as inputs, and with Motion 
Durations, Submotion Durations, Collision Branch Lower Limits, Sequence 
Lower Limits, also as inputs, the Optimize Unwrapped Schedule 88 portion of the 
Computerized Model 64 will determine whether There is A Feasible Schedule? 
5 90. As shown, there is an additional input: Thermal Forming Process Duration 
"N", which includes the Target (time), Limits and Lock Status. 

The operator may, for example, decide that a defect is occurring because 
there isn't enough "reheat" time and input a proposed new reheat time. The 
operator could also input more than one new Thermal Forming Process Durations 
10 N1,N2,..., during an offline evaluation of the process. In either of these modes 
the Event Angles for the entire schedule would be available and these could all be 
inputted by the operator or downloaded from the control for the machine. 

If no schedule is feasible, the Computerized Model will Reject The Inputs 
92. If a schedule is feasible, the Computerized Control will Output The Thermal 
15 Forming Process Durations 89. Such an output might, for example, be a printout 
for each duration, of the target duration, an indication of whether or not its target 
duration was locked, and the actual duration located in a window extending 
between the high and low limits for the duration. Should there be a solution the 
Wrap Optimized Event Times Into Event Angles 84 portion of the Computerized 
20 Model converts the Event Times to Event Angles and proceeds to Print Event 
Angles And New Machine Cycle Time 94. 

Figure 19 illustrates the use of the computerized model for complete 
schedule optimization (Schedule Optimization). Machine Cycle Time, Event 
Times, Motion Durations, Submotion Durations, Thermal Forming Process 
25 Durations, Collision Branch Durations and Sequence Branch Durations which 
represent target values, are possible inputs to Optimize The Unwrapped Schedule 
96. In addition, a number of limits are also inputs: 1. Min/Max Motion Duration 
"N", 2. Min/Max Thermal Forming Process Duration "N", 3. Min/Max Collision 
Branch "N", and 4. Min/Max Sequence Branch "N". The Min/Max Motion 
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Duration "N" relates to servomotor driven displacements which can be selectively 
varied. Given these inputs, The Optimize Unwrapped Schedule finds and 
optimized schedule if a feasible schedule exists. In the event that the query There 
Is A Feasible Schedule ? 98 is answered in the negative, the operator will be 
5 advised to Loosen Limits 100 so that the operator will try to find a solution by 
modifying the limits. In the event that the query There Is A Feasible Schedule? 
98 is answered in the affirmative, the control may Set Collision/Sequence 
Branches To Max, Lock All Other Durations and Again Optimize The 
Unwrapped Schedule 101. This will maximize these branches to further reduce 

10 the rate of collision of missequencing. The computer model will then Wrap 
Event Times Into Event Angles 102, Print The Event Angles And The New 
Machine Cycle Time 104 and Output Optimized Durations VS. Limits 106. The 
operator accordingly has the ability to manipulate the unwrapped schedule to the 
fullest. He can start with an existing job file which traditionally would have the 

15 cycle time, event angles, and the servo-motion branch durations and work to 
define an optimized schedule. Alternatively, he could enter the Thermal Forming 
Process Durations and convert them to the Event Times (a screen, not shown 
could display all of this information to facilitate his analysis). 

The Computerized Model can, if There Is A Feasible Schedule 107 

20 (Figure 20), determine if There Is An Active Constraint(s) That Restricts Further 
Improvement? 108 and will Output The Active Constraints) (including the 
direction to move for improvement) 110. For example, the computerized model 
may show that the constraint that is preventing optimization is blow mold cooling 
time. This then enables the operator to address this specific problem to increase 

25 the flow of coolant through or at the molds. If there is no solution the operator is 
advised to Loosen Limits 100. 

Figure 21 illustrates the use of this technology to optimize the wear on a 
mechanism operated by a servomotor (Wear Optimization). Here the 
Computerized Model 64 is used to optimize an unwrapped schedule and when 
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There Is A Solution 107, the next step is for the computer model to Optimize The 
Unwrapped Schedule Locking All Variables Except Servo Motion Durations And 
Setting Target Servo Motion Durations At Large Value. 1 12. The next step is for 
the computer model to print The Optimized Duration For Servo Motor "N:" and 

5 to Deliver The Optimized Duration For Servo Motor "N" To Servo Motor "N" 
Controller 114, which will then Route Duration Of Servo Motor "N" From Servo 
"N" Controller To Servo "N" Amplifier Drive Card 1 16 which will then Change 
To Optimized Duration In Amplifier Digital Signal Processor 118. The 
Amplifier Digital Signal Processor could, for example, scale a normalized motion 

10 profile for the mechanism to be driven to accommodate any duration of motion. 
In this environment, an ideal motor to adjust in this fashion is a servomotor which 
has a normalized motion profile that could be scaled from a minimum duration to 
a maximum duration. While the preferred embodiment of profiled activator is a 
servo motor, other electronic motors, such as a stepping motor, could be used. 

15 The disclosed control can be used with a glass forming machine either 

directly as a part of the machine control or indirectly as a control of a machine 
which is virtually operated for evaluation purposes. 
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